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PREFACE 


The low degree and order terms in the spherical harmonic model of 
the tidal potential have been observed through the perturbations which 
are induced on near-earth satellite orbital motions. This recovery, 
which is the most complete dynamic model ever obtained has been 
achieved through evaluating tracking observations on 17 dl . ffer J2, 
mostly laser, satellites. A new improved GEM-T1 geopotential model, 
complete to degree and order 36, was estimated simultaneously with the 
66 adjusted tidal coefficients. The gravitational and tidal models were 
developed using the J2000 Reference System with the adopted nutations of 
Wahr and the precession model of Lieske. The tidal recovery was made in 
the presence of an extended oceanographic model containing over 600 long 
wave] ength coefficients from 32 major and minor tides. Since solid 
earth tides have identical perturbing frequencies as the ocean tide. , 
the solid earth tidal model of Wahr was used as a basis for the recovery 
of the ocean tidal terms. This provided a complete description of the 
combined tidal potential sensed by these well tracked satellites. This 
tidal model (for all 32 adjusted and unadjusted tides) has then been 
used to calculate the secular change in the moon's mean motion due to 
tidal dissipation and the tidal braking of the earth’s rotation. The 
secular change in the moon's mean motion due to tidal dissipation is 
found to be -25.27 ± 0.61 arcsec cy . Our estimate of the luna 

acceleration agrees well with that observed from Lunar La«|r Ranging 
techniques, which most recently found -24.9 ± 1.2 arcsec cy (Newhall 
et al 1986). The corresponding tidal braking of the earth s 
is -5.98 ± .22x10 -22 rad sec -2 . If the non-tidal braking of the earth 
due to the observed secular change in the earth's second zonal harmonic 
(Yoder et al , 1983) is considered, modern satellite techniques yie a 
total value of the secular change in the earth’s rotation rate of - . 9 
± 0. 36x1 0 -22 rad sec -2 . 
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1.0 INTRODUCTION 


The gravitational attraction of the sun and moon deform the earth, 
its oceans and its atmosphere. These deformations, known as tides, 

provide a rich spectrum of effects. Among those which are readily 
observable, the tides cause temporal variations in the external gravita- 
tional attraction the earth exerts on near-by orbiting objects and cause 
complex changes in the rate of the earth's rotation. 

Any imperfect tidal response of the earth results in a tidal 
bulge which is not exactly aligned in the direction to the disturbing 
third body. This phase lag in the tidal response produces a torque 
which causes a transfer of angular momentum in the earth/moon/sun 
system. Since the ocean tides usually have much larger phase lags than 
the solid earth tides, they play a dominating role in this exchange of 
angular momentum. This transfer of angular momentum manifests itself as 
secular changes in the orbital elements of the disturbing bodies and in 
the braking of the earth’s rotation rate. 

The subject of tidal dissipation is a classical problem which is 
discussed at length in Lambeck (1980). Tidal dissipation has been a 
subject of intense interest in the geophysical, astronomical and oceano- 
graphical communities, because the redistribution of the earth's mass 
due to the tides is the recognized explanation for the observed accel- 
eration in the earth's rotation and acceleration of the lunar orbit. 

The earth/moon and earth/near-earth satellite systems have strong 
physical similarites. Man-made earth orbiting satellites have certain 
observational advantages over the moon. They have better known, and in 
most cases, laboratory measured physical characteristics. Man-made 
orbiting objects are generally easier to observe with great precision by 
earth-based ranging systems. The diversity of objects available for 
tracking also helps in the separation of contending forcing effects. 
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Therefore, precise measurements to artificial satellites can play an 
important role in understanding the evolution of the earth/ moon/sun 
system. 


One disadvantage of resolving the tides from their effects on 
artificial satellites lies in the short observation history available on 
these orbits and the resulting inherent difficulty of separating secular 
from long periodic effects. However, the near earth satellite record is 
lengthening in time, permitting more robust tidal solutions. Centimeter 
level laser systems now routinely perform both lunar and near - earth 
satellite tracking. Direct comparisons and inferences can be made of 
the common satellite responses--either that of the moon or an artificial 
sat el 1 it e — - to the external potential of the earth, especially those 
resulting from tidal forces. 

Theoretical results associated with the evolution of the earth/ 
moon system have not previously been in close agreement with those 
obtained observationally . This disagreement has caused some to introduce 
other physical mechanisms, such as a temporal change in the gravita- 
tional constant, to bring theory and observation into balance (Van 
Flandern, 1975). Recoveries of tidal parameters from artificial satel- 
lites give important new insight into the earth's tidal response and the 
role these deformations play in the evolution of the lunar orbit. 

The early work of Kaula (1964, 1969) first addressed the effects 
of solid earth tides on near-earth satellite orbits. Lambeck et al 
(1974) complicated the picture by showing the importance of ocean 
tides. Additionally, Lambeck (1977, 1980) has discussed a wide range of 
geophysical and astronomical consequences of these tides, such as their 
importance in the evolution of the earth/moon system, and braking the 
rotation of the earth. The purpose of this paper is to evaluate, in 
these contexts, the consequences of a new satellite based solution of 
earth and ocean tide parameters. 
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Most of the early analyses for the recovery of tidal parameters 
from orbital perturbations were investigations of the long period tidal 
effects evidenced in the evolution of the orbit's mean Kepler elements. 
In most cases, the inclination and ascending node of the orbit were 
analyzed. These two out-of-plane components of the orbit have strong 
tidal signals, and these components are less sensitive to errors in 
other force models, such as the non-conservative effects. Tidal values 
from satellite analyses have been appearing in the literature since the 
late 1960's. Some of the major satellites used in our solution were 
previously analyzed in this manner; Douglas et al (197*0 and 

Felsentreger et al (1976) studied Geos-1 and Geos-2, Cazenave and 
Daillet (1981) performed an analysis of Starlette, Goad and Douglas 
(1977) studied an Oscar Navy navigation satellite and Goad (1977) 
evaluated Geos-3. 

The introduction of very precise laser ranging systems and an 
accompanying strong global tracking network made the combined recovery 
of tidal and geopotential terms possible from a direct analysis of 
satellite tracking data. As will be shown, the tides and the geopoten- 
tial are closely related so that their effects can both be represented 
in the usual spherical harmonic form (Larabeck, 1980). Herein, the 
gravitational model is treated as the constant part of a time varying 
potential, with time dependent effects solved for at tidal frequencies. 
Williamson and Marsh (1985) broke important new ground in this direction 
with their tidal solution from Starlette and Christodoulidis et al 
(1986) followed with a comparable analysis of Lageos. 


In this paper we report on a new dynamical tidal solution which 
has been obtained at NASA/Goddard Space Flight Center as part of an 
effort to develop an improved geopotential model (Marsh et al, 1987) in 
preparation for the 1991 launch of Topex/Poseidon (Born et al , 1986). 
In order to yield an accurate time-invariant gravity model, a complete 
separation of gravitational and tidal potential is necessary. This was 
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accomplished through improved modeling, including the adoption of 
special tidal models which were more complete than those ordinarily 
available from oceanographic sources and the simultaneous least squares 
estimation of numerous geopotential and tidal coefficients. 
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2.0 MODEL 


A natural form for expressing the gravitational field for orbital 
calculations is in terms of spherical harmonic coefficients, 


V 



“ 4 a 4 

m-0 (f£) P H«l <Sln * ,[C ta 0 ° S “ X * S lm Sln miil 


( 1 ) 


where GM is the earth's gravitional constant (including the atmosphere), 
a e is the earth's mean equatorial radius, is the associated Legendre 
function of degree 4 and order m, and r,<|> and \ are the distance from 
the center of mass, latitude and longitude (e.g. Kaula, 1966). A number 
of the potential coefficients C„ and S„ are considered to be time 
invariant unknowns and are allowed to adjust based upon observed orbital 
motion. 


It is shown (Appendix 1) that in a system where the disturbing 
bodies (sun and moon) are referred to the ecliptic, and the perturbed 
body (the artificial satellite) is referred to the equator, the second 
degree tide potential is given by: 
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for the solid earth and by 
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for the oceans. This spherical harmonic form provides a linear response 
at each constituent f in the tide- generating potential. The tidal and 
gravitational potentials combine to perturb the orbit. As can be seen 
in the above equations, the second degree tidal potential is equivalent 
to that due to gravitation with tidally driven, temporally varying 
second degree harmonics. Aliasing of the gravitational solution can 
result if there i 3 a poor temporal sampling of these tidal effects 
and/or if the tides are poorly modeled. The principal earth and ocean 
tides, their Doodson numbers and appropriate Darwinian correspondancies 
are shown in Table 1 . 

The principal perturbations due to solid earth and ocean tides 
have identical frequencies on satellite orbits. Consequently, earth and 
ocean tidal parameters cannot be recovered separately from the analysis 
of near-earth tracking data. In our solution, an adjustment is made to 
a set of ocean tidal parameters, which corrects the observed temporal 
exterior potential of the earth relative to an adopted solid earth tide 
model and relative to a partially unadjusted base ocean tidal model. By 
design within the solution, a sufficient number of independent tidal 
terms have been adjusted at their dominant frequencies. This solution 
and its interpretation form the basis of this paper. 


2.1 A PRIORI SOLID EARTH AND OCEAN TIDE MODELS 

For the solid earth tides, the frequency dependent model developed 
by Wahr (1979) has been adopted and held unadjusted within the solution 
(Table 2). All of the remaining solid earth tides are modeled with a 
frequency independent Love number of k2 = «30 with zero phase. Their 
tidal potential is evaluated in the time domain from lunar and solar 
ephemerides . 
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Schwiderski (e.g., 1983, 1981, 1980) and Parke (1982) have 

analyzed global data acquired by tidal stations and have estimated 
oceanographic models for the major ocean tides. These models have been 
spectrally decomposed using quadrature under the assumption that the 
tidal amplitude (A) and phases (iJj) are known at all points on the globe 
(where A 3 0 for land areas). The A^cos^j, and Aj.sini^ are computed on a 
uniform dense grid of points on the earth’s surface to solve for 
the a dm,f’ b dm,f’ c Hm,f and d dm,f used t0 ^present each tide in 
harmonic form: 


A f cos$i f 
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(4) 


For each tide constituent (f), the values used for orbit computations 
are obtained as: 


C dm,f COS l *dm,f “ 2 ^ a dm,f " d dm,f^ 

C dm,f Sin C dm,f " 2 ^ b Slm,f + C dm,f^ 

(5) 

C dm,f COS C dra,f 2 " ^ a JLm , f " d dm,f^ 

C dm,f Sin C dm,f \ ~ [ °dm,f + b dm,f ] 
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The expansion is taken out to some maximum order, itmax. The ocean 
tidal phases in the oceanographic models can be related to that of eq. 3 
by: 


e - - -5 + tt 6 

f f 2 om 


( 6 ) 


where 6 om is the Kronecker delta (where m=0 denotes the long period, m-1 
the diurnal and m=2 the semi-diurnal tides) and 


X * 


0 

IT 


positive Doodson coefficients 
negative Doodson coefficients 


(7) 


Note that for our major tides, x =1T only for K1. is used in this form 
in Appendix 1 . 

There are many minor ocean tides which are not developed directly 
from oceanographic data due to the computational burden required for 
precise solutions of the Laplace Tidal Equations. Because of the 
coupling of the satellite motion and the disturbing body, many of these 
tides have significant effects on satellite orbits and are a source of 
aliasing in the recovery of satellite derived ocean tidal models if not 
properly accounted for. This is especially true of the numerous 

sideband constituents which typically have an amplitude of 10$ of the 
major tide at a nearby frequency. Total neglect of these sideband tides 
significantly degrades the accuracy of the adjustment for the terms 
defining the major tides. To calculate models for many minor ocean 
tides which were otherwise unavailable, an algorithm was developed and 
tested. Figure 1 diagrams the predictive method which was developed for 
the computation of any tide constituent and its errors given its 
frequency. The algorithm uses linear admittances to predict any tide in 
the long period, semi-diurnal and diurnal bands (Christodoulidis et al , 
1986a). 
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A sensitivity analysis was performed using the linear orbit 
perturbation theory of Appendix 1 to isolate those tides which were 
important for the satellite solution, and ocean-based tidal models were 
developed for those found to be satellite sensitive. The long wavelength 
coefficients of the ocean tides were adjusted simultaneously with the 
geopotential if sufficient sensitivity existed within the satellite 
observations. Otherwise, the ocean tidal terms were held fixed at their 
oceanographically determined values. Note that a change in the solid 
earth model is accomodated by a comparable change in the solved for 
ocean tides of the same frequency. This is further discussed in 
Section 3. Table 3 summarizes the complete set of ocean tidal terms 
which were modeled in the GEM-T1 solution and indicates which terms (if 
any) were adjusted within each tidal frequency. 

Due to attenuation with altitude, most near-earth satellites are 
insensitive to ocean tidal terms above degree 6 and are most sensitive 
to the long period perturbations which arise exclusively from the 
prograde terms in the harmonic expansion. R. Eanes of the University of 
Texas has shown (private communications) that there are a few ocean 
tidal terms which we have neglected. These neglected terms are not of 
the same order as the dominant terms and resemble m-daily gravity 
effects of a f ew cm amplitude in orbital perturbations. As m-daily 
effects, these neglected tidal terms can cause an aliasing error in the 
geopotential recovery, but they are quite distinct from the tidal 
parameters adjusted herein. 


2.2 THE GEM-T1 SOLUTION DESIGN 

The recent GEM-T1 solution (Marsh et al., 1987, in press) using 
significantly more and better tracking data was a complete re-iteration 
of Goddard Space Flight Center's general gravitational model solutions. 
The availability of a super computer (the Cyber 205) and the 
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vectorization of our major software tools greatly impacted the timeframe 
required to complete this work. GEM-T1 was carefully designed to be a 
general gravitational field (i.e., not a "tailored field" for any 
specific orbit). GEM-T1 was solved in the J2000 Reference System (with 
IAU 1980/Wahr nutation and modern reference constants being adopted). 
This reference system was generally based on the Merit Standards 
(Melbourne, et al., 1983). A summary of the reference parameters 
adopted within GEM-T1 can be found in Table 4. 

The GEM-T1 solution made use of a large subset of the precise 
third generation laser tracking data (with tracking precision better 
than 5cm) taken on 7 geodetic satellites. Strong Tranet Doppler data 
sets from Seasat and Oscar 14 were also used. Interestingly, each of 
the geodetic satellites which have been studied earlier for tides have 
now been reanalyzed as part of this simultaneous multi-satellite 
gravitational and tidal solution. In each case, improved tracking data 
sets (i.e., better precision and global distribution) have been utilized 
in the formation of GEM-T1 . Furthermore, a simultaneous adjustment of 
tidal and geopotential terms has yielded more complete and accurate 
dynamic tidal results. The GEM-T1 gravitational model is complete in 
spherical harmonics to degree and order 36; 66 tidal parameters 
representing the longest wavelength coefficients of 12 major tides were 
simultaneously adjusted in the solution. Approximately 550 other ocean 
tidal terms were used to model long wavelength ocean tidal variability 
for 32 major and minor tides. 

In all, 17 satellites were included in the GEM-T1 solution. Of 
the 17, weaker data sets from six low inclination objects were used 
specifically to help resolve the zonal harmonics of the field. The 
strong data sets were those acquired by recent laser and Doppler 
systems. Table 6 and Figure 2 describe the orbital characteristics of 
the satellites used in the formation of GEM-T1 . 
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The summary of the data in Table 5 reveals that precise laser 
tracking played a dominant role in defining the GEM-T1 geopotential and 
tidal models. Of particular importance, two laser satellites, Lageos 
and Starlette, made strong contributions to both the tidal and 
geopotential field recoveries. These are completely passive orbiting 
objects which were exclusively designed to serve as space-based laser 
targets. Both satellites are extremely dense spheres (area to mass 
ratios of .00069 and .00096 m 2 kg -1 respectively), covered by laser 
corner cubes, and are in orbits designed to minimize non-conservative 
forcing effects. Lageos orbits at nearly an earth radius above the 
earth and senses only the longest wavelength gravitational and tidal 
fields. Starlette, orbiting at a much lower altitude of -1000 km, 
experiences a rich spectrum of tidal and gravitational perturbations and 
is highly complementary to Lageos for the separation of long and short 
wavelength gravitational and tidal terms. Figures 3»1 and 3.2 show the 
complimentary tidal perturbation spectra sensed in the inclination and 
ascending node perturbations of Starlette and Lageos. Both of these 
satellites are tracked on a high priority basis by a global network of 
laser tracking stations and have extensive observation sets which have 
been acquired under NASA's Crustal Dynamics Project activities (Coates 
et al, 1985), Project MERIT, and the WEGENER Campaign. 

A summary of the principal periods of the solid earth and ocean 
tides on the major satellites utilized in the solution can be found in 
Figure 4. These tides perturb the satellite orbits with dominating long 
period effects. On the earth's surface, the largest tides are found in 
either the diurnal or semi-diurnal band. The precession of a satellite's 
ascending node and not the earth's rotation is responsible for the wide 
range of these satellite periodicities. 
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3.0 THE GEM-T1 OCEAN TIDAL SOLUTION 


Table 3 has summarized the ocean tidal terms which were selected 
for adjustment in the GEM-T1 solution (that is, adjusted simultaneously 
with the gravitational field). Tables 7.1 through 7.4 present the 
recovered ocean tidal spherical harmonic coefficients respectively for 
each degree. The values shown for the coefficient and phase uncertain- 
ties were taken from a covariance matrix which was calibrated for the 
geopotential coefficient uncertainties (see Marsh et al, 1986). We 
believe this calibration also applies for the tides. This calibration 
gave realistic (one sigma) uncertainties for the estimated tidal and 
geopotential terms. Because the ocean tidal terms are inseparable frcm 
those of the solid earth at the same frequencies, the uncertainty 
specified for a given tidal term does not reflect changes which would 
arise from modifications to the fixed solid earth tides which are 
modeled. Rather, these uncertainties represent the error in the sum of 
the solid earth and ocean tidal models at a given frequency mapped into 
the spherical harmonic decomposition of the ocean tide. These uncer- 
tainties represent commission errors (in a strict sense) in modeling the 
external potential of the earth at tidal frequencies when using the 
complete ocean tidal description available in GEM-T1 plus the Wahr solid 
earth tidal models. Neglect of an adjustment for the very long period 
tides (e.g. the 18 year tide) is a potential source of aliasing which is 
not considered in the error estimates which are presented. 

Note that the GEM-T1 recovered tide coefficients do not require 
correction for sideband effects as was necessary in Cazenave & Daillet 
(1981) and in Williamson and Marsh (1985). The errors associated with 
the sideband constituents have been minimized because the sideband terms 
are explicitly modeled in our solution and the satellite data spans a 
significant number of years. 


13 


preceding page blank not filmed 


PAGfc_ / WAN* 


In the current solutions, the Wahr Earth Model is adopted to 
describe the solid earth's tidal response. This model was developed 
assuming a uniformly rotating earth comprised of an elastic inner core 
and mantle and a liquid core. Dehant (1986a, b,c) has estimated the 
effects of anelasticy on different components of the earth's tidal 
response using models developed by Zschau (1985) and others. She has 
shown that the neglect of anelasticity may be the cause of some of the 
discrepancy in the observed versus predicted gravimetric factors at 
tidal frequencies. Wahr and Bergen (in press) have also addressed the 
question of the effect of anelasticity on the earth tides (and in the 
tidally induced variation in the earth's rotation rate) and have 
estimated upper and lower bounds for these effects. They have noted 
that observational data of sufficient precision is lacking to permit a 
meaningful estimate of anelastic parameters directly, even though recent 
recoveries of nutation amplitudes by space based technologies like VLBI 
put us on the threshold of doing so. Wahr and Bergen (ibid) show that 
for the diurnal and semi-diurnal tidal frequencies, the lower and upper 
bound for the anelastic effect on k^ is typically .002 and .015 
respectively. These assessments agree well with those developed by 
Dehant. Using the midpoint of these upper and lower bounds shown by 
Wahr and Bergen for the anelastic contribution at specific tidal 
frequencies and a 1$ change in i <2 for all other tides in our model, we 
have performed a sensitivity study by re-estimating all 66 ocean tide 
coefficients in the presence of this changed solid earth tide model. 
The adjusted 2nd degree ocean tidal coefficients which result generally 
change in phase by approximately 3 degrees and in amplitude by less than 
10$. However, the resulting external potential described by the 
complete ocean and solid earth model remains nearly constant (within one 
third of the standard deviations given). This test of the stability of 
the model's external potential compared the calculated values of the 
lunar acceleration and the braking of the earth's rotation with the 
values obtained with the original Wahr model. 
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The GEM-T1 satellite determined values are also compared to those 
estimated from the oeanographic tidal solutions of Schwiderski and Parke 
(Tables 7.1 to 7.4). Generally, the comparison of satellite and oceano- 
graphic tidal values is quite good. The variation seen between the two 
oceanographic tidal solutions is often larger than the uncertainty of 
the dynamic satellite solution. We believe these results confirm that 
the adopted Wahr solid earth tidal model is quite reasonable and that 
the solution for long wavelength ocean tidal terms is realistic when 
compared to those observed oceanographically. These satellite determined 
ocean tidal coefficients can provide effective constraints within future 
oceanographic-based model developments. 

Of special interest are the second degree values for the M 2 tide 
where the amplitudes show large disagreement between oceanographic 
models. For this term, which dominates the lunar evolution, the GEM-T1 
satellite result favors the value obtained by Parke. For the 0 1 second 
degree values, a large disagreement is also present, only this time the 
GEM-T1 satellite results favor the Schwiderski model. 

The large uncertainty in the oceanographic models (as evidenced 
by the level of disagreement they exhibit) prevents an accurate 
separation of solid earth and ocean tidal contributions in the observed 
total effects seen by orbiting objects. As explained in Williamson and 
Marsh (1985), significant improvements are also required in the modeling 
of atmospheric tides, atmospheric pressure loading over both the oceans 
and the continents, the radiational potential, and other such subtle 
effects. This lack of strong earth-based modeling is the reason that 
satellite tidal analyses cannot yet be expected to contribute 
significantly to assessing tidal dissipation in the solid earth. Zschau 
(1986) comes to the same compelling conclusion in his attempt to 
interpret the M 2 solid tide dissipation from combined satellite-ocean 
tide solutions. 
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Table 8 compares 2nd and 4th degree ocean tidal coefficients 
obtained by other satellite studies with the values from the GEM-T1 
solution. Our satellite values seem to be consistent with most 
previously reported satellite results, especially the M2 analyses. For 
M 2 , both the second and fourth degree terms are in good agreement. Note 
the Cazenave and Daillet (1981) values for 0-| shown were taken from 
Williamson and Marsh (1985), who corrected them to refer to the Wahr 
nutations from the original values using the Woolard nutations. 

The GEM-T1 ocean tidal model is a significant improvement over 
earlier studies. It alone is based upon the simultaneous adjustment of 
geopotential and tidal terms, it uses the Wahr nutations, it has modeled 
550 presently unrecoverable sideband and higher degree ocean tidal 
terms, and has been derived from complete tracking information on a 
multiplicity of satellites. Most, if not all, of these factors are 
significant limitations in the earlier satellite ocean tidal 
investigations. Because of the nature of present model improvements, 
the systematic errors in the earlier investigations are probably 
comparable to or larger than any error estimates found in previous 
analyses. We have therefore not tabulated the comparitive error 
estimates from these earlier investigations. 
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4.0 SECULAR PERTURBATIONS IN THE MOON/SUN ORBIT 


The observed secular acceleration in the moon's orbit is dominated 
by tidal dissipation within the earth's oceans. The secular variation 
of the Kepler elements of the moon and apparent sun are calculated using 
the tidal potential in the Lagrange Planetary Equations, evaluating 
those terms of zero frequency as detailed in Appendix 2. The ecliptic 
Kepler rates take the form: 


da* 

- K* 

+ 

sin 

+ 

dt 
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2m, f 

2m, f 

de* 

= K* 
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sin 
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dt 
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= K* 
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sin 
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dt 

i 

2m , f 
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( 8 ) 


The * stands for either the moon or the sun (there is also a small 

change in the apparent solar orbit); m = 0 for long period tides (LP), 

m=1 for diurnal tides (D) and m=2 for semi-diurnal tides (SD). These 

+ + 

labels are used in summary tables. and are the second degree 

tidal amplitudes and phases obtained from the GEM-T1 tidal solution or 
the background unadjusted tidal models, and K*, K*, and K* are 

3 0 1 

calculated from analytical theory as shown in Appendix 2. 

To obtain the third body accelerations from the GEM-T1 tidal 
model, it is necessary to consider the contributions of all of the 
adjusted and unadjusted second degree terms. The Wahr solid earth tidal 
model has no anelasticity and thus there is no phase lag in the earth's 
lunisolar tidal response and no resulting torque on the lunar orbit. 
The complete second degree ocean tide model includes the unadjusted 
tidal coefficients derived from admittances (Table 9) which have been 
included in all of our near-earth orbital calculations as well as those 
terms adjusted from the satellite data (as given in Table 7.1). 
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The secular changes in the Kepler elements of the moon and 
apparent solar orbit based on the GEM-T1 model due to tidal energy 
dissipation in the earth's oceans are summarized in Table 10. The effect 
of the earth's ocean tides on the lunar orbit is reasonably well 
determined and can be compared with lunar evolution observations. 
However, only the change in the apparent solar semi-major axis is 
measured significantly in the sun/earth case. The uncertainties in the 
Kepler-element rates which are given have been obtained from a formal 
RSS propagation of the tidal coefficent errors (as given in Table 7.1) 
combined with a 10$ estimated uncertainty for all unadjusted terms shown 
in Table 9. 

Using the simple relationship derived from Kepler's Law, the 
secular perturbation of the lunar mean motion is obtained from the semi- 
major axis rate as 

dn ~3 n* da 

dt 2a* dt 


where the computation of n*, the mean motion, is given by 


n* 



and 


M* 


M f 1 



( 10 ) 


(ID 


Details of the contribution of each tidal constituent and its associated 
uncertainty to the secular change of the moon's mean motion is given in 
Table 11. The total acceleration from the GEM-T1 ocean tidal model is 
the sum of each of these contributions and, as shown in Table 12, 


18 


produces a value of -25.27 ±0.6 arcsec/cy 2 for the secular change in 
the moon's mean motion. This value is in excellent agreement with latest 
value obtained observationally by Lunar Laser Ranging: (LLR) (Newhall et 
ai., 1986 in press) of -24.9 ± 1 arcsec/cy 2 . The n m observed from LLR 
is in harmony with the energy dissipation calculated from the earth's 
ocean tidal model for the earth/moon system. This confirms that energy 
dissipation due to the ocean tides is the dominant source of the lunar 
orbit's secular acceleration and no other physical mechanisms are 
required at the current .6 arcsec cy~ 2 level of precision. 

Table 11 also reveals that the long period tides have secular 
effects on the lunar mean motion. Our formulation allows us to evaluate 
these effects: the coefficients describing the tides are equatorially 

referenced whereas the satellite orbital frequencies of the tidal 
disturbing potential are best described and implemented herein as linear 
combinations of the ecliptic rates. Previous derivations have missed 
this phenomenon because of their expansion within an equatorial refer- 
ence system. In this equatorial system, the correspondence between the 
tidal perturbation rates and the orbit element rates do not exactly 
cancel, whereas the more correct derivation (given in Appendix 2) find 
them doing so, giving rise to secular effects. 

In our calculations we have ignored any contribution that arises 
from the dissipation of energy due to tides raised on the moon by the 
earth. While M2 is not of consequence due to the permanent bulge it 
causes in the moon's shape, contributions from 0 1 and N 2 may need to be 
considered. Lambeck (1975) estimates that these terms have small 
effects, being no more than .16 m/cy in — • for the moon. This effect is 
of the order of the accuracy we are now achieving and requires further 
investigation. 
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5.0 TIDAL BRAKING OF THE EARTH’S ROTATION 


The ocean tides retard the axial rotation of the earth. Conserva- 
tion of angular momentum in the earth/moon/sun system requires that the 
secular change in the rotational velocity of the earth due to the 
lunisolar tidal torques is (Bursa, 1986): 


dn 

V — m 
1 dt 


+ A 


de 

n 

2 dt 


m 


+ A 


di 
n 

3 dt 


m 


+ B 


dn 

c 

1 dtT 


+ B 


de 

c 

2 dt" 


+ B 


di 

c 

3 dt" 


0 dC 
C dt 


( 12 ) 


where a , e , i are the ecliptic Keplerian elements of the 
m m m 

moon; a , e , i are the equatorial Keplerian elements of the apparent 

3 3 3 

sun; the A terms are for the orbital angular momentum contributions of 
the Earth/Moon system; the B terms are for the orbital contributions of 
the Earth/Moon system about the sun; and the C term represents the 
effects due to changes in the Earth’s polar moment of inertia C. 

The formulae for the A and B terms, which are taken from Lambeck 
(1980), have been derived from conservation of angular momentum in the 
direction of the Earth’s mean spin axis. In Lambeck’ s derivation, he 
uses the inclination of the lunar orbit on the equator averaged over one 
period of the lunar ascending node on the ecliptic, I m . Thus 

cos I * cos i cos i 
m s m 
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The requisite formulae are then, where M is the mass of the Earth, M s is 
the mass of the Sun, and M m is the mass of the Moon: 


m 


M M 


1 n a cos i 

m 


mm s 

M+M 

m 


C 


, m 

A, » •= — cos 1 
1 3n m 

m 


A_ = F e cos i 
2mm m 


A_ - F sin i 
3 m m 


for the lunar terms, and 
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for the solar terms. 
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Bursa (1986) has evaluated the effects due to changes in the 
Earth's polar moment of inertia due to the secular time variations in 
the second zonal harmonic: 


§ - I M&2 tht kg m2/cy 

dt 3 at 

where Ma 2 is 2. 4296x1 O 38 kg m 2 . The corresponding C 
and Q is 7.2921 15x1 0" 5 rad/ sec. Yoder et al (1983) 


( 21 ) 

is 8.0378 x10 37 kg m 2 
have estimated: 


UJ ^ — Q 

— - = (-2.8±0.6) x 10 /cy 
dt 


from Lageos tracking data, and have attributed it to the effects of 
post-glacial rebound. Thus, the non- tidal acceleration of the earth is 


* _ fi dC 

n NT = ' C dt 


-22 2 

(1.29±.28) x 10 rad/sec 


Each second degree tide in the GEM-T1 expansion, including both 
adjusted and modeled values, has been used to calculate its contribution 
to braking the earth's rotation. These rates with their associated 
uncertainties are given in Table 11. A summary of the long period, 
diurnal and semi-diurnal contributions to fl and that from the non-tidal 
secular change in J 2 is presented in Table 12. Our estimate of the 
secular tidal acceleration of the earth's rotation (n T ) is -5.98 ± 
0.22x1 0 -22 rad sec -2 . Taking into account the non-tidal effect of the 
change in the earth's oblateness (J 2 ), our best estimate of the secular 
braking of the earth's rotation is -4. 69^± 0. 36x1 0* 22 rad see' 2 . In 
terms of the change in length of day, 0 T corresponds to 2.25 ± .08 
msec cy - ^ and 0 corresponds to 1.76 ± .14 msec cy 

Lambeck (1980) summarizes the astronomical results which have 
been obtained for constraining the values of Q and n m and for the 


23 


determination of h itself. Although there appears to be a wide varia- 

tion of values, Lambeck found h = -5.48x10 -22 rad s“ 2 to be his 
consensus value from astronomical observations. Our estimate of this 
rate, which was satellite derived using both a new ocean tidal model and 
the observed non-tidal change in the earth's 2nd zonal harmonic, agrees 
reasonably well with his interpretation of the astronomical constraints. 
Clearly, this new satellite determination lacks many of the problems 
associated with the determination of the time and place of ancient 
eclipses, and represents a significant advancement in determining the 
tidal braking of the earth’s axial rotation. 

The derived value for Si discussed herein, is solely based on the 
long term secular effect. While the tidal contribution, n ? , is complete 
and not expected to vary over the three thousand years of recent history 
(Lambeck, 1980), the non-tidal contributions may not be complete. In 
our derivation of 0 NT we have based the estimate on the effects which 
arise from J 2 . Yoder et al (1983) suggest that other smaller 
contributions may come, for example, from J 3 , although present orbital 
histories cannot accurately resolve them. 

Recently, a re-analysis of the eclipse reoords by Stephenson and 
Morrison (1984) suggests that there is evidence of a non-uniform 
decrease in the earth's rotation rate over the last three thousand 
years. However, the physical basis for any such change remains 
unexplained. In the case of Stephenson and Morrison, the evidence for a 
changing non-tidal acceleration is based upon an apparent discrepancy 
between the historical eclipse results prior to 1620 and the analysis of 
the later telescopic data subsequent to 1620. From their telescopic 
data, they derive changes in the length of day of 1.4 msec cy -1 versus 
2.4 msec cy 1 for the earlier eclipse data. The artificial satellite 
results for fl T , the tidal effect alone, yield 2.25 ± .1 msec cy -1 , which 
is a value quite consistent with that Stephenson and Morrison found from 
their early eclipse records. Our value of 1.76 ± .14 for Si is not really 
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too different from their 1.4 given the accuracies of the astronomic 
analyses. Their telescopic results depart from those found herein using 
modern satellite laser ranging techniques by about -.4, and the early 
eclipse results depart by about .6. A simple explanation for these 
disagreements may be that the astronomic data is of insufficient 
accuracy; also error sources such as inconsistent reference frames may 
be present. Given the lack of supporting error analyses, their 
conclusion of non-uniformity in the decrease of the earth's rotation 
rate is questionable. 

From another perspective, Lambeck (1980, eqn. 10.5.2) estimated 
that the tidal deceleration of the earth's rotation and the secular 
change in the lunar mean motion are related by 

« T =* (51 ±4) n m 

This is simply derived from the separate estimates of n m and If we 

use this equation with our estimate of n m , we find 6^ = -6.27x1 0 - ^ 

rad/sec , which is in good agreement (one sigma) with our solution. Our 

• « 

estimates of n m and 0^. provide the corresponding relationship 
fi T “ ( 49 ±3 ) ^ 

indicating that we find the Lambeck formula slightly overestimates the 
effect of rijj, on ftp. 

These results from the present artificial satellite analyses 

• • 

represent the Q and n m of an epoch around 1980. These results can be 
directly compared to the excellent contemporary observation history of 
the earth's rotational behaviour available from the space technologies 
of Satellite and Lunar Laser Ranging and Very Long Baseline Inter- 
ferometry to assess the contemporary role tides play in this complex 
retardation of the earth's rotation rate. 
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6.0 SUMMARY 


Modern laser tracking technologies allow us to monitor the orbital 
evolution of near-earth satellites to unprecedented levels of accuracy. 
Recently acquired observation data sets from a total of seventeen mostly 
laser satellites have been used to develop detailed models of the forces 
which explain these observed orbital effects. Both the time invariant 
and the temporally varying tidal gravitiational potential fields of 
the earth have been accurately resolved from our analysis of direct 
satellite tracking observations. This model, Goddard Earth Model (GEM) 
T1 , is the most complete gravitational and dynamic tidal solution ever 
produced exclusively from direct satellite tracking data. 


The GEM-T1 solution simultaneously recovered a complete gravita 
tional field to degree and order 36 and 66 tidal coefficients repre- 
senting twelve major tides. The solution made use of modern geodetic 
constants, such as the J2000 Reference System complete with the Wahr 
nutation and solid earth tidal models. A more complete model of the 
background ocean tides was used, which had been developed to reduce 
aliasing in the recovery. These unadjusted tidal terms included over 
550 coefficents representing the longest wavelength constituents of 
32 major and minor tides. The minor ocean tidal models were estimated 


using a predictive technique based on the observed major tides and a 
postulated linearity of the admittances. The tidal perturbations which 
result from the separate solid earth, atmosphere and oceans give rise to 
identical orbital frequencies which cannot be separately resolved frcm 
an analysis of orbital motion. However, the total tidal effect is 
accurately sensed in these well tracked orbits and the tidal 
perturbations are known to be strong. 


The complete GEM-T1 tidal model containing both the adjusted and 
unadjusted terms has been used to estimate the effects of the earth’s 
tides on the earth/moon/sun system. The secular change in the lunar mean 
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motion estimated from this dynamic tidal model, -25.27±.6l arcsec cy’ 2 , 
agrees very well with the value of -24.9 ± 1. arcsec cy~ 2 recovered 
from a decade of Lunar Laser Ranging by Newhall et al , 1986. The tidal 
braking of the rate of the earth’s rotation due to conservation of 
angular momentum in the earth/ moon/ sun system yields a value of 
^•9 ± .22 x 10 rad sec , corresonding to 2.25 ± 0.08 nisec/cy in 

length of day. Thus we have derived a set of tides from tracking of 
artificial satellites which seems to consistantly predict the observed 
evolution of the lunar orbit. 

Taking into account the non-tidal braking due to the Lageos 

derived secular rate of change in the earth's second zonal harmonic 
(Yoder et al , 1983), artificial satellite techniques provide a total 
value for the secular change in the earth's rotation rate of 
-4.69 ± .36 x 10 22 rad sec 2 , which corresponds to 1.76 ± .14 msec/cy 
in length of day. This value is consistent with recent astronomic 
studies of eclipse records such as that of Stephenson and Morrison 
(1984), who obtained 2.4 msec/cy for the data prior to 1620 and 1.4 
msec/cy for the telescopic data since 1620. Because there may exist 
other secular effects which we have not yet accounted for, and because 
of the significantly poorer accuracies acheivable frcm the astrcmical 

data, this agreement must be regarded as quite good. Also, the 

1.76 msec/cy does not confirm any inference of a non-uniform change in 
the earth's rotation rate; rather it tends to indicate that the 
difference of 1 msec/cy in the separate estimates of Stephenson and 

Morrison may reflect the separate errors of their determinations. 
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TABLE 1 


Darwinian 

Symbol 

Doodson' s 

Argument 

Number 

Period 

(hr) 

Description 

m 2 

255.555 

1 2. 42 

Principal lunar semidiurnal 

s 2 

273.555 

12.00 

Principal solar semidiurnal 

n 2 

245.655 

12.66 

Larger lunar elliptic 
semidi urnal 

k 2 

275.555 

1 1 .97 

Lunar/Solar semidiurnal 

l 2 

265.455 

1 2. 19 

Smaller lunar elliptic 

K 1 

165.555 

2 3.93 

Lunar/Solar diurnal 

°1 

145.555 

25.82 

Principal lunar diurnal 

P 1 

163.555 

24.07 

Principal solar diurnal 

M f 

075.555 

1 3. 66d 

Lunar fortnightly 

M m 

065.455 

27.55d 

Lunar monthly 

S sa 

057.555 

1 88. 62d 

Solar semi-annual 
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TABLE 2. WAHR LOVE NUMBERS MODELED IN GEM-T1 


Band 

Tidal Line (f) 

k 2 



±t_ 

Long Period 

All 

.299 

Diurnal 

145555 (01) 

.298 


163555 (PI) 

.287 


1 65545 

.259 


165555 (K1 ) 

.256 


165565 

.253 


166554 (PSI) 

.466 

Semi-Diurnal 

All 

.302 


*A11 phases are equal to zero. 
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Table 3 


OCEAN TIDE MODELING 
FOR 

GRAVITY RECOVERY 


• LONG PERIOD TIDES • 


Ooodson 

No. 

Darwin 

Name 

Modeled 

Adjusted 

056.554 

Sa 

deg. 2-*6 

deg. 2 

057.555 

Sgo 

prograde 

deg. 2 

058.554 


only 

none 

065.455 

075.555 

075.565 

^ ID 

M, 



deg. 2 
deg. 2 
none 



• DIURNAL • 



135.655 

Q t deg. 2-*6 

none 

145.545 

prograde 

none 

145.555 

Oy and 

deg. 2, 3, 

155.455 

retrograde 

none 

155.655 

M, 


none 

162.556 

Pit 


none 

163.555 

P, 


deg. 2, 3, 

164.556 

Si 


none 

165.545 


none 

165.555 

*i 


deg. 2, 3, 

165.565 


none 

166.554 



none 

167.555 



none 

175.455 

Si 


none 

185.555 

00, J 


none 


36 



• SEMI-DIURNAL • 


Doodson 

Darwin 




NO. 

Nome 

Modeled 

Adjusted 

245.655 

N 2 

deg. 2- *6 

deg. 2, 3, 4, 5 

255.545 


progrode 

none 

255.555 

m 2 

and 

deg. 2, 3, 4, 5 

265.455 

1 2 

retrogrode 

none 

271.557 




none 

272.556 

h 



deg. 2, 3, 4, 5 

273.555 

5 2 



deg. 2, 3, 4, 5 

274.554 

R2 



none 

275.555 

K2 



deg. 2, 3, 4, 5 

285.455 




none 

295.555 




none 
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Table 4 


CONSTANTS FOR NORMAL 
EQUATION GENERATION FOR 
GEM-T1 


APR I OR I GR AVITY MODEL 


GEM-L2' FOR 
PGS-I33I’ FOR 
PGS-S4’ FOR 
GEM- I OB’ FOR 


LAGEOS 

STARLETTE 

SEASAT 

ALL OTHER SATELLITES 


PRIME C) INDICATES FIELDS RESOLVED 
HOLDING C,S 2 , = 0 

EARTH TIDFS 


WAHR WITH k 2 = 3.0° PHASE 

nominal 


CHRISTODOULIDIS ET AL MODEL WITH 
600 INDIVIDUAL TERMS, PREDICTED FROM 
OCEANOGRAPHIC MODELS USING 
ADMITTANCES 

TIDAL DEFORMATIONS 

h 2 AND A 2 AT MERIT VALUES 

GM, j, ,_LZL 

GM = 396600 436 km’/sec 2 
o. = 6376 137 km 
1/f = 296 257 


POLAR MOTION AND AI-IJTI 

ZERO MEAN SET OF POLAR MOTIONS 
BASED ON 6 YEARS OF SLR DETERMINATIONS 


STATION c oordinates (unadjusted) 

250 STATIONS TRANSFORMED TO LAGEOS 
SL-6 AND ROTATED INTO NEW POLAR 
MOTION SYSTEM 

THIRD BODY EFFECTS 

PLANETS MERCURY THROUGH NEPTUNE 


Z AXIS 

INSTANTANEOUS SPIN-AXIS 


J2000 


RELATIVITY 

NONE 
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Table 5 


DATA UTILIZED HI PRELIMINARY 
TOPEI GRA VITY MODEL: 1986 


NUMBER OF 

SATELLITE DATA TYPE TOMALMAI&1CIS- 


LAGEOS LASER 

STARLETTE 

GEOS-l 

GEOS-2 

GEOS-3 

BE-C 

SEASAT 

Dl-C 

Dl-D 

PEOLE * 


56 

46 

48 

28 

36 

39 

14 

4 

6 

6 


SUB-TOTAL - LASER 


285 


SEASAT DOPPLER J5 

OSCAR-14 * 

SUB-TOTAL - DOPPLER 28 


GEOS-l CAMERA 

GEOS-2 

ANNA 

TELSTAR 

BE-C 

COURIER IB 

VANGUARD-2RB 

VANGUARD-2 

Dl-C 

Dl-D 

PEOLE * 


43 

46 

30 

30 

50 

20 

10 

10 

10 

10 

9 

6 


SUB-TOTAL - CAMERA 


273 


TOTAL 


560 * 


•PEOLE arcs contained both optical and laser data. 


NUMBER OF 
OBSERVATIONS 

144527 

57356 

71287 

26613 

42407 

64240 

14923 

7455 

11487 

4113 

444.408 


136042 

63098 

201.140 


60750 

61403 

4463 

3962 

7501 

1739 

2476 

686 

1299 

2712 

6111 

38 

153.140 


798.688 
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Table 6 


SATELLITE ORBITAL CHARACTERISTICS 


SATELLITE 

SATELLITE 

NAME 

ID NO 

ANNA-18 

620601 

BE-B 

640841 

BE-C 

650321 

COURIER- IB 

600131 

Dl-C 

6701 1 1 

Dl-D 

670141 

CE0S-1 

650891 

CEOS-2 

680021 

CE0S-3 

750271 

LACE0S 

760391 

OSCAR 

670921 

PE0LE 

701091 

SEASAT 

780641 

STARLETTE 

750101 

TELESTAR- 1 

620291 

VANGUARD-2RB 

590012 

VANGUARD-2 

59001 1 


SEMI-MAJOR 
AXIS 


7501 

7354 

7507 

7469 

7341 

7622 

8075 

7711. 

7226 

12273 

7440 

7006 

7170 

7331 

9669 

6496 

8296 


ECC 

.0082 

•0135 

.0257 

.0161 

0532 

.0848 

0719 

.0330 

0008 

0038 

.0029 

.0164 

0021 

0204 

2429 

1832 

1641 


INCL 
(DEG ) 

50 12 
79 69 
41 19 
28 31 
39 97 
39 46 
59 39 
105 79 
114 98 
109 85 
89 27 
1501 
108 02 
49.80 
44 79 
32 92 
32 89 


DATA* 

TYPE 


* D e Doppler 
L t Laser 
0 1 Optical 
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Table 6 

COMPARISON OF OCEAN TIDAL COEFFICIENTS 

FROM 

SATELLITE SOLUTIONS 


M 


m 


M 


f 




n 2 

m 2 


s 


• • 

2 


K 


2 


Amp 

Phase 

Amp 

C * 

U 2m 

£ * 

C 2m 

C* 

4m 

0 36 

274 10 


1.99 

132 5 


1 .80 

245.35 


1.70 

265 4 


2.69 

318 53 

1 .83 

2.83 

327 

(1.43)* 

2 85 

344.3 

1 .44 

2 96 

354. 

1 .7 

0.81 

296.83 

0.35 

0.73 

253. 

(0 63)* 

1 77 

1 68.3 

0.74 

2.61 

328.50 

2.62 

2 63 

31 17 

3.40 

0 70 

334 0 

0 17 

0.59 

346. 

(0 37)* 

0.76 

322 9 

0.24 

3.26 

320.93 

0.93 

3 45 

319. 

(1.01)* 

2.57 

319 2 

1 .70 

3.5 

338 

1 .03 

3.1 

333. 

1 03 

3 4 

325 

1 .0 

3.2 

331 . 

0.9 

2.9 

332. 

10 

0.80 

301 .93 

0.41 

0 77 

302. 

(0 21)* 

0 78 

32 4 

0 42 

0.31 

302 14 

0.19 

0 16 

238 

(Oil)* 

0.83 

23.4 

0.53 


Phase 

Z\ 

4m 

Source 

(ref) 


GEM-T1 (paper) 
STARLETTE (B) 


GEM-T1 
STARLETTE (B) 

283.38 
(276.3)* 
304 5 
350 0 

GEM-T1 
LAGEOS (C) 

STARLETTE (B) 
STARLETTE ♦GEOS 3 (D) 

234 91 
(258.3)* 
264.9 

GEM-T1 
LAGEOS (C) 
STARLETTE (B) 

254 42 

244.4 

GEM-T1 

STARLETTE (B) 

137 29 
(103.0)* 
192.7 

GEM- T 1 
LAGEOS (C) 
STARLETTE (B) 

127.41 
(124.7)* 
1 18.5 
79. 

78 

124. 

113 

115. 

GEM-T1 
LAGEOS (C) 

STARLETTE (B) 
STARLETTE (E) 
STARLETTE ♦GEOS 3 (D) 
OSCAR, STARLETTE, 
GEOS 3 (F) 

GEOS 3^0SCAR (G) 

GEOS 3 (H) 

86 53 
(141 .8)* 
94.8 

GEM-T1 
LAGEOS (C) 
STARLETTE (B) 

75.96 

(103.5)* 

96.0 

GEM-T1 
LAGEOS (C) 
STARLETTE (B) 


••includes atmospheric tide 
* held unadiusted at Schwiderski value 

(B) Williamson and Marsh, 1985 (F) Fclsentreqer et al., 1979 

(C) Christodoulidis et al., 1985 (G) Goad and Douglas, 1978A 

(D) Cazenave and Daillet, 1981 (H) Goad and Douglas, 1978B 

(E) Cazenave and Daillet, 1981 
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Table 9 


VALUES FOR MODELED 2ND DEGREE TIDES 
CALCULATED FROM LINEAR 
ADMITTANCES 


TIDE (f) 


+ 

C o 

2m, f 

AMPLITUDE ( cm) 



PHASE ( deg) 


LOWC PERIOD (m = 0) 

058.554 

075.565 


0.03 

0.35 


231.63 

261.15 


DIURNAL Cm 

= D 



135.655 

«1 

0 53 

31 3.70 

145.545 

Off 

0.47 

31 3.85 

155.455 

"if 

0.06 

31 4.03 

155.655 

"l 

0.17 

314.04 

162.556 

*1 

0.05 

31 4.22 

164.556 

Si 

0.02 

31 4.26 

165.545 

K 1f 

0.05 

31 4.28 

165.565 

*1S 

0.36 

31 4.28 

1 66.554 

Vi 

0 02 

31 4.30 

167.555 

®i 

0.04 

31 4.32 

1 75.455 

Ji 

0.13 

314.60 

185.555 

00) 

0.06 

315.08 


SEMI-DIURNAL (m = 2) 


255.545 

M 2$ 

0.1 1 

265.455 

*-2 

0 07 

271.557 


0 002 

274.554 

*2 

0.01 

285.455 


0.01 

295.555 


0.002 


316.82 

315.86 

314.93 

314.57 

312.22 

308.07 
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Table 10 


SECULAR CHANGES IN THE ECLIPTIC KEPLER 
ELEMENTS OF THE MOON AND THE SUN 
DUE TO TIDES 


dt 


de« 

dt 


.die 

dt 


da G 

dt 


de Q 

dt 



3.73 ± 0.09 


(1.83 ± 1.10) x 10 


(-6.65 ± 0.53) x 10 


(1.43 ± 0.25) x 10 


(1.3 ± 9.4) x 10 19 


(-5 ± 51) x 10 


m cy 


-1 


yr 


-i 


deg yr 


mey 


yr 


deg yr 
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Table 1 1 


SECULAR CHANGE IN THE MEAN MOTION OF THE 
MOON (n. in ARCSEC/CY 2 ) AND IN THE 
ROTATIONAL VELOCITY OF THE EARTH 
(0, in 10" 22 RAD/SEC 2 ) DUE TO TIDES 


TIDES 


CJn 

— “<•: 

Q 

LONG PERIOD TIDES 

056 554 S a 

0.00 

0.00 

000 

057 555 S sa 

- 0.00 

0 00 

0 04 

056 554 

- 0.00 

0.00 

- 0.00 

065 455 M m 

-0 03 

0 31 

-0 00 

075 555 M f 

-0 56 

0 16 

-0 12 

075 565 

-0 10 

0 06 

- 0.01 

TOTAL LP -TIDES 

-0 69 

±0 36 

-0 17 


DIURNAL TIDES 
135655 0 , 

-0 16 

0 04 

- 0.03 

145 545 0 1f 

-0 10 

0 02 

-0 01 

145 555 o , 

-2 92 

0 25 

-0 65 

155 455 M 1f 

0.00 

0.00 

0.00 

155 655 M, 

0 01 

0.00 

0.00 

162 556 H \ 

-0 00 

0.00 

- 0.00 

163 555 Pi 

- 0.00 

000 

-0 12 

164 556 S , 

000 

0.00 

0 00 

165 545 K 1f 



- 0.00 

165 555 K,moon 




165 565 K,sun 




165 565 Kis 


— 

0 01 

166 554 Y 1S 

0.00 

0.00 

0.00 

167 555 Oi 

000 

0.00 

0 00 

175 455 J , 

0 01 

0 00 

0.00 

165 555 00 , 

-0 00 

0 00 

-0 00 

TOTAL D-TIDES 

-3 16 

±0 25 

-0 60 




0.00 
0 03 
0.00 
0 12 
0.07 
002 


±0 14 


0 01 
0 01 
0.07 
0.00 
0.00 
0.00 
0 10 
0.00 
0.00 


0.00 
0 00 
0 00 
000 
0 00 

± 0.12 
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♦ 


TIDES 

Ho: 

SEMI-DIURNAL TIDES 

245 655 N 2 

-1 43 

255 545 M 2S 

0 02 

255 555 Mj 

-20 00 

265 455 L 2 

0 0 1 

271 557 

-0.00 

272 556 T 2 

-0 00 

273 555 S 2 

-0 00 

274 554 R 2 

0 00 

275 555 K 2 m°°n 


275 555 K 2 sun 


265 455 

-0 00 

295 555 

0.00 

TOTAL SD -TIDES 

-2 1 40 





Q_ 

Sfl 

0 16 

-0 2 1 

0 04 

0 01 

0.00 

0.00 

0 40 

-4 45 

0 09 

0.00 

0.00 

000 

0.00 

- 0.00 

0.00 

0 00 

- 0.00 

0 05 

0.00 

-0 35 

0 04 

000 

0 00 

0.00 


0.00 

- 0.00 

0.00 

0.00 

0.00 

000 

±0 43 

-5 oi 

±0 12 
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Table 12. 

SUMMARY OF THE SECULAR TIDAL ACCELERATION 
OF THE EARTH Q AND THE ACCELERATION OF THE 
MOON n FROM THE GEM-TI OCEAN TIDAL MODEL 


6(10 22 rad/sec ~ 2 ) 

6 T Lp -0.17 ±0.14 

G T D -0 60 ±0.12 
sd -5 01 ±0.12 


n (arcsec cy~ 2 ) 

n T Lp ”0 69 ± 0.36 

n T D -3. 18 ±0.25 
n T SD -21.40 ±0.43 


6 t -5.98 ± 0.22 n T -25.27 ±0.61 

Q hj *1.29 ±0.28 


Q -4.69 ± 0.36 


n -25.27 ±0.61 
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BACKGROUND TIDES MODEL DFVEl.OPMFNT 



X 


tidol heights (A), 
phases (y) and 
errors given by 
Schwiderski for 
major tide 
constituents 
l 

COMPUTE 




I 

| = known admittances for f j , f ^ . fj frequencies 

(computed from Schwiderski 
1 tides) 

| = predicted admittance for fp’s 

^ compute 

heights (A) , phases (y) & errors 
for fp’s 

* weighted least squares 


FIGURE t 
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FIGURE 2. ORBITAL CHARACTERISTICS OF THE GEM-T1 SATELLITES 





056 5545 

057 5555 

065.4555 

075 5555 
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163.5555 

££££'£91 

m 
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in 

so 

in 
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272.5565 
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275.5555 

SATELLITE 

B 

s,. 


« f 

D 

P 1 

<1 

D 

m 2 

D 

s 2 

m 

LAGEOS 

365 

183 

27.6 

13.7 

13.6 

221 

1050 

9.20 

14.0 

159 

280 

524 

STARLETTE 

— 

365 

163 

27.6 

13.7 

1 1.9 

60.8 

91.0 

7.61 

10.5 

33.1 

36.4 

45.5 

GEOS- 1 

365 

183 

27.6 

13 7 

12.6 

65.4 

160 

6.20 

1 1.7 

48.3 

55.7 

80.2 

GEOS- 2 

365 

183 

27.6 

13.7 

14.4 

629 

257 

9.83 

15.3 

2250 

436 

129 

GEOS- 3 

365 

183 

27.6 

13.7 

15.2 

482 

132 

m 

17.2 

145 

104 

66.2 

BE-B 

365 

183 

27.6 

13.7 

13.1 

118 

332 

8.66 

12.6 

70.2 

RQ 

166 

BE-C 

365 

183 

27.6 

13.7 

1 1.8 

57.9 

84.8 

7.51 

10.3 

31.5 

34.4 

42.4 

SEAS AT 

365 

183 

27.6 

13.7 

14.8 

7130 

178 

10.2 

16.1 

331 

174 

m 

TELSTAR- 1 

365 

183 

27.6 

13.7 

12.8 

93.9 

193 

8.34 

12.0 

53.9 

63.2 

96.7 

ANNA 

365 

183 

27.6 

13.7 

12.0 

64.4 

99.4 

7.71 

Qg 

35.3 

39.1 

49.7 

OSCAR 

365 

183 

27.6 

13.7 

13.6 

180 

11700 

9.12 


119 

177 

5830 


Figure 4 
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APPENDIX A 


Mathematical Formulation 
of the 

Farth and Ocean Tidal Potentials 


The purpose of this appendix is to provide the mathematical 
formulation of the solid earth and ocean tidal potentials appropriate 
to the use of osculating ephemerides. The relationship between the 
Kaula type indices and the equivalent Doodson numbering scheme is 
provided Additionally, the linear orbit perturbations for the Kepler 
elements are given 


The formulation of the solid earth and ocean tidal potentials 
for use in Cowell type programs implies that the satellite accelerations 
are computed as the gradient of the potential: 


r r ((DAr) = 


dU 

dr 




3U 

ad) 


rydi,v) = 


r cos (J) 


3U 

a \ 


> 


In our software, 6E0DYN, the accelerations are evaluated 
in the inertial system aligned with the instantaneous earth-fixed 
system These accelerations are then rotated to the inertial system 
of integration via transformations accounting for precession, 
nutation, polar motion and earth rotation 


( 1 . 1 ) 


These tidal potentials are formulated in terms of the 
ecliptic elements of the disturbing bodies and the equatorial elements 
of the satellite. The development we followed is quite similar to Goad 
(1977), wherein the surface harmonics due to the disturbing body 
are rotated from the equatorial to the ecliptic system This 
formulation thus can use the planetary ephemerides (currently 
DE-200) directly, and the tidal perturbations are computed with the 
identical precise ephemeris used for other perturbations Using 
this form means that there are no terms which are dependant on 


A-l 


both the sun and moon: these interaction effects are automatically 
generated by the numerical intergration process. It should also be 
noted that the distinction between mean end osculating is not currently 
significant (R J. Eanes, 1981, private communication). 

T he Ex pr essions De scribin g The Solid Earth And Ocean Tide Potentials 
The following initial definitions are required: 


(J), A, r 



are the latitude, east longitude, and radial 
distance of the point of evaluation, 

is the Greenwich sidereal hour angle. 


Fjq (cos O') 


are the associated Legendre functions 


a *, e*, i*, 
0*. M* 


are the Keplerian elements of the disturbing 
body referred to the ecliptic. 


€ m , 'm 
^m, w m, 


are the Keplerian elements of the moon referred 
to the ecliptic. 


a s 


As 



are the Keplerian elements of the sun referred 
to the ecliptic. 




m 


is the gravitational constant times the mass of 
the moon 




is the gravitational constant times the mass of 
the sun 


is the gravitational constant times the mass of the earth 
is the average radius of the earth 


A-2 


is the equivalent of the Doodson constant, 


G d 

2 

f 


indicates summation over all tide constituents 
(f s) in the expression of the tide generating 
potential 


A is the equivalent of the Doodson coefficient^ 

k is the second degree love number and <^ 2 ^ 

2 '* is its phase, 


k are the load deformation coefficients, 

9 

p 0 is the average density of the oceanic water, 



are the amplitude and phase in the spherical 
harmonic expansion of the ocean tides 
specified by Jl, q, ± , and the tide 
constituent f 


It can be shown that, in a system where the disturbing bodies 
(sun and moon) are referred to the ecliptic whereas the perturbed 
body is referred to the equator, the second degree tide potential 
is given by 

V(d) . A. r )= 2 k 2 ,f G d ("jr) P 2 m ( n $ ) cos cx ^ 


for the solid earth, and by: 

U((t.A,r)=2 I 4nGRPo (1^ ) c: , (sin®) 

f M,± 1 V 2JI+ 1 / r Jlq 


C08O< PqJ (1 3) 


for the oceans, where the angular arguments are respectively 


<x* = ( + ) 


(2-2h)oo* + (2-2h+j)M*+k Q' 


] +m V 


m 


* n-m n / 2 ♦ 


( 14 ) 


and 
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moon 


sun 


= ( + ) (2-2h)co* + (2-2h+j)M*+kO*J ♦ m0 4 ± q A 


Tf - m ^ / + £ 

'2 HA 


( 1 . 5 ) 


Note that there is a difference between ± which belongs to the 
ocean tide expansion and (±) which is part of the definition of the 
tidal constituent f The tidal constituent f corresponds to 
specification of m, k, h, j, (±), 8nd *. The equivalent of Doodson s 
constant is given by 


G = 4 




3 


4 fl 3 
8 m 


( 16 ) 


The equivalent of the Doodson coefficient is given by 




( 17 ) 


where 


A is the Kronecker delta 

om 


d = A ok + 7 1-k Wk 


(1 8 ) 


C =6-1 4j T 


1 \ ^ r(±) 

2mk 


1 

(2-k) i 
(2+k)i 


(-) 

(♦) 


( 19 ) 


t, <± . ) = (co 3 n m(i)k Ginn (±)k_m F (±) (co 3 2 o 

2mk 2mk 


(1 10 ) 


F^ (x) = - d — [" 2 -m, x _ |) 2 ' , ' rn 1 

2mk W <j x 2(±)k L X ( x W J 


(1.11) 


P 

l - -y- and £ is the obliquity of the ecliptic. 
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The inclination function F 2kh (i*) 8nd eccentricity function G 2h j(e*) 
can be found in Kaula (1966, pp. 34-39). 

The formulation of the potential in (1.2) and (1.3) is compatible 
with Doodson's development (Doodson, 1921) Doodson, however, 
elected to adopt Brown's lunar theory for the description of the 
position of the disturbing bodies, and thus his formulation is in 
terms of the mean elements of the sun and moon. In the present 
development, osculating elements are used 

The Relationship To Doodson A g ument Numberin g 

It is often convenient to express the angular arguments in the 
same form as the Doodson solution. This requires the following 
osculating angular variables which correspond to Doodson's mean 
variables: 

l - local osculating lunar hour angle, measured from lower transit 
of the moon past the local meridian 

l- <x - s - ti 

where o< is the right ascension of the point of evaluation and s is 
the moon's osculating longitude given by 

s = O m * * M m 

p = moon's osculating longitude of perigee 
P ~ + CO m 

N"= the negative osculating longitude of the moon's 
ascending node 

N'= -O m 

h z the sun s osculating longitude 
h z Q s + co s + M s 

p ( = the osculating longitude of the sun's perigee 
p, = Q s * co s , and 

N^z the negative osculating longitude of the sun's 
ascending node 

n;= 
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(1 12 ) 

(113) 

(114) 

(1 15) 
(1 16) 
(1 17) 

(1.18) 



The angular argument of the tidal potential used by Doodson, 
allowing for the use of osculating variables, is: 

A p = p,C ♦ p 2 s * p 3 h ♦ P 4 p PsN' + P 6 P, + P 7 n; 

Where p, 's are integers which are typically in the range of -4 to 4. 
As p, is always positive, the angular argument was expressed by 
Doodson as an "argument number". 

P, (P 2 *5) (P 3 + 5)* ( p 4 ^ 5) ( p 5 5 ) ( pe + 5 ) ( P? + 5 ) 

Note that p 7 is not present in Doodson's development because the 
mean longitude of the sun's ascending node is by definition zero 
on the ecliptic 

Using (1 19) and (1 I 2), the angular argument takes the form 

A p = P,o< + (p 2 - p,)s * p 3 h+ p 4 p + p 5 N' ♦ p 6 p,+ p 7 N, - p,n 

Ommission of the phase <^ 2 ,f in equation (1.4) produces the 
angular argument of the tide generating potential which is 
equivalent to that expressed by equation (12 0 except for a 
possible offset. Allowing for the following definitions of the 
astronomical longitudes of the disturbing bodies 

L* = co * + M* ♦ Q * 

P* = 00 * + Q* 

the angular argument in ( 1 4) takes the form 

& 2 ,\ = me* ♦U-mj (?) £(2-2h+j) L* - jP* -(2-2h-k) Q* 

Note that L* is either s or h, P* is either p or p,, and Q* is either 
-N'or -N, Also, the entire argument is either for the moon or for 
the sun, as mixed terms are not present in this formulation 
(Doodson has a small number of mixed terms which are present 
due to his choice of independent variables ) 

By equating the coefficients of the independent variables in 
(1 23) and (121), the following correspondences are obtained 
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...(1.19) 


...( 1 . 20 ) 


( 121 ) 


(1 22 ) 


(1 23) 



a) Lunar Case 



Pi = m 

p 5 = (-) (2-2h-k) 

(1.24) 

P 2 -P 1 = (*) (2-2h+j) 

Ps = 0 


= 0 

p4 = (*)j 

P 7 = 0 


b) Solar Case 

Pi = m 

O 

II 

in 

QCL 

(1 25) 

E 

II 

CCL 

II 

CN 

CEL 

Pt = ( s )j 


P 3 = (?) (2-2h+ j) 
P 4 = 0 

p 7 = (±)(2-2h-k) 


The inverse correspondences are: 
a) Lunar Case 


Lei k = p, - P 2 - P 4 + Ps 
For k > 0: 


(1.26) 

m = Pi 

J - P 4 


h = (P 2 + P 4 - Pi ♦ 2)/2 
k = k 

sign = (+) 


For k < 0 

m =p, 

J = - P 4 


h = (p^ - P 2 " P 4 + 2)/2 
k = -k 

sign = (-) 
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b) Solar Case 


Let k _ p 7 - Pe - P 3 


(1.27) 


For k > 0 

m = Pi J - Ps 

h = (p 3 + p G *2)/2 sign = (+) 

k = k 


For k < 0 

™ = Pi j = - P* 

h = (-p 3 - p 6 + 2)/2 sign = (-) 


k = -k 


Note that when either k or m is zero, there is a twofold 
symmetry in that two different choices of indices generate the 
same frequencies When both m and k are zero, a fourfold 
symmetry exists This is so whether the Doodson notation or the 
kaula type indices are used 

The Doodson coefficients indicated by ( 1 .7) must be summed 
over the symmetric cases to compare with Doodson's tables, in 
these cases the values of the separate coefficients are identical 
Additionally, when the angular argument is independent of the 
elements of the disturbing body (h=1, j=0, k=0), each disturbing 
body has a twofold symmetry and the coefficients must also be 
summed over the disturbing bodies (i e , for the permanent tide, 

K, , and K 2 ) 


ORfQfNAL PASE IS 
OF POOR QUALITY 
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Perturbations In The Kepler Elements 


The analytical formulations for the tidal perturbations of the 
Kepler elements of a satellite are derived below. As can be seen 
from inspection of equations (12), (13), (1.4) and (15), the 
solid earth potential is a special case of the ocean tide potential 
Thus the following derivations are for the ocean potential, and 
the solid earth case will be derived therefrom. 


.et 




i!L GR P . o ( iTii \ r + 


p 


2q J 


^ 5 q,f °Sq,f + d A, and 


(1 28 ) 
( 1 . 29 ) 


]J - the gravitational constant times the mass of 
the earth 


The tide potential 
be rewritten as: 


U,(UU,r)= 2 K,;, 

M,± q ' J 


given in ( 1.2) for the oceans can then 
£ 

7T+T P fiq (sjn ^ C0S ( Q £qJ ± qA) 


( 1 . 30 ) 
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Letting 


Nq ,f = K SqJ C0S Q SqJ + K Pq,f C ° S 


(131) 


B 


*qj = Sm CT lVf + ^ 


♦ K" i sing 


$q,f 


particular 8,q term of the potential can be written as 


u >,.r 


A Pq j cos qX ♦B Jq ^ sin qX 


(1 32) 


which is of the same form as the potential given by Keula (1966 


p 31, eq (3,53)) Carrying on in 


the same fashion as Kaula did, (1 30) 


takes the f orm (see ibid eq 3,70 & 3,7 1 ): 


U #Qi =-^_ 2 F (i) 2 Gp ps (e)S Jqps f (go, M,0, 0 g ) 

a$+i p =o vqp 3 = - <x) 


(1 33) 


where a, e, i, co , M, O are the Kepler elements of the perturbed 
body, the functions F pqp and G pps ore given in Kaula (1966, 
pp 34-38), and 


S *qps,f ( ^ M ' 0 ' 1 V-- 


A pq.f 

'B pq J 


P -q even 


P-q odd 


cos x pqps + 


® j>q-f 

■# 4 > ^ J p-q odd 


P -q even 


sin x ( 1 34) 


pqps 


with 


X 


pqps 


2p ) go + ( Jl ~ 2p + s)M + q (Q - 0g) 


(1 35) 


It can be easily shown that S pqps j may also be written as 


SpqpsJ (W,M,O,0g) = | K^j 


COS 

sin 


P -q even 
P -q odd 


( X pqps ± *^pq,f) (136) 
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The time rate of change of the argument 


*^pqps,f X 9qps i ® tq.i 

is dominated by the equatorial secular rates in CJ, Q , and M of 
the satellite, by the rotation of the earth represented by 0$, 
and by the ecliptic secular rates of the disturbing body in 
<^>*,Q* # 8nd M*. Hence YpqpsJ is approximately constant over 
the period associated with Sj, qps j. The linear perturbations in 
the orbit elements can be constructed (Kaula, 1966, pp 37 - 41) 
for each term desired by analytically integrating the Lagrangian 
equations of motion for that term 


a _ llP fl 2 F flqp G pps (A- 2p+s) 

A °wJ ' WR 




Y Pqps J 


(1 38) 


Ae , F »gp G »Ps( 1-e ? )^[(1- 

pqps,f M 


e 2 ) 2 (A- 2p+3 MJl - 2p ) 




s 


pqps,f 


(1 39) 


A 


sqps,f - 


F pqp G pps [(A- 2p ) cos i-q] 
na p + 3 (1-e 2 )2 sin i Y~ qps f 


^pqpsj 


(1 40) 


AQ ( = pR 

pqps,f 


(^ F pqp^. ) Gp ps Sj qps j 


na 


9-3 (i- e 2 )2 sin i ¥; qps j 


+ bQ 


(141 


A<o qps ,f = pR J [( 1 -e 2 )2 e -'F Sqp ( dG »% e ) - cot i ( 1 -e 2 n( aF «%, ) G Sps ] 

na fl + 3 YpqpsJ 


G sqps ,f 
* bud 

(1.42) 


AM 


pqps,f 


= HR .[-(l-e 2 )e-' ( aG »-^, e ) + 2(JI*0G, ps - 3nG, ps (J -2c, )(f J qpj ] 


hQ P 3 ¥ pqpsj 


(1 43) 


• Fp qp Sjj qps j + j ri 
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where n is the mean motion of perturbed body. 



SjqpsJ is the integral of with respect to Us argument: 

and <5ft, Aco>, AM are the indirect effects due to the earth's 
oblateness From Balmino (1974): 


sin 

•cos 


ji-q even 


£-q odd 


T ^ qps J 




r n 

i 



SQ 


20-e 2 ) 2 


4 cos l 

A CJ 

= 3nC 20 (|7 

i 


e 


(i - e 2 ) 2 | 

i 

|0-e 2 ) 

1 “5 C03 1 

AM 


1 

3 


1 . 

1 - 3 cos 2 l 



? 3. 

4 (l-e 2 )^ 




...(1.44) 

(1.45) 


(1 46) 


+ 


-si n i 


4 si n 2 1 


si n 2 i 


with the earth's second degree harmonic C 20 ~ - 1 .082628* 1 0 3 
and 


E 1 ACpqps J ( t ) 
l - 1 pqps ,f ^ ) 


(1 47) 


Integrating the inclination, we get 

I = jjR 


Fjqp G Sps [01- 2p) cos i - q] 


na ,+3 (1-e 2 )* sin i (j J qps j) 


2 ^Jqps.f 


(1 48) 
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and for the eccentricity 


r _ ^ qp (I - ® 2 ) 2 [( 1 -e 2 ) 2 (Jl- 2p+3) -(Jl- 2p)j - 

" ) 2 S '- J 

Substituting (1.48) and ( 1 49) to ( 1 46) and by letting 

q , _ 4UGR 2 p 0 f 1+kp n 

M ps 'f p \ 2 jl + V Uj!qps ^ 


and 


^ JlQPS A 

4rrGR p 0 ^ 

f l+k * ' 

) ^jqpsJ 


\2Jt ♦ 1> 

where 




^ jiqps A 

= y r- 
t S<U 

COS 

sin 

p-q even 

p-q odd "ft pqps J 

^ oqps J 

- y r± 
t Sq,f 

sin 

-cos 

p-q even 

P ~ q ° dd TiqpsJ 


we have 

= 6ttG Pq G2q ( _R ) ( 1^1 \ F<iqp g hps q 
<l-e*P». a k 2 * ♦ I ' ft* ,) 2 #,PSJ 

V U Jlqpsj / 

J/ 

• q + 4 cos l ( 1 -e 2 ) 2 ( jl -2p + s)- 5 cos i ( Jl -2p) 


b 


6nGp 0 C 20 / r \ 5+4 f 1 + k, 


co - TT. 

(1 -e 2 ) 2 a 


(-*) ( 1 ** ) fi g G * P s q 

0 v 2Ji *] J \ 2 ^ ps '< 


(if flqpsj ) 

• |^5 cos 2 i ( Jl -2p) - 5q cos l ♦ 2(1-5 cos 2 l )( 1 -e 2 )' 2 ( Jl -2p+s)-2( Jl -2p) | 


1 


. (1 49) 

(1 50) 

(151) 

..(1.52) 

(1.53) 

(1 54) 
...(1.55) 
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= 


9UGp 0 C 2 Q 
(l-e 2 ) 2 a 


r R_^* 4 / 'l*kp 

1 a J ♦! 


) 


Fjiqp Gj ps p 
/ • ± \ 2 jqpsJ 

IT sqpsj / 


|5 cos 2 i a-2p>- 2q cos i -(Jl-2p) + (1-3 cos 2 i) ( 1 -e 2 P Ul-2p+s) 


From (1 38) through (1 56) we can write the final form of the 
perturbation in the Kepler elements due to the oceanic tidal 
potential 


Aa x = 
«qps \ 


• ± 

n Y" 


r ^ s+2 f 

2 jt ♦ i t pqp “Pps 


cf) ( 


) F M p G w JI - 2p * 5)c w.f 


sqpsj 


_ 4ti Gp 0 ( 1 -e 2 )' 2 / R \ 9 * 2 / l*kj| \ 

^ e pqps,f • + a ^ 2il + 1 pq p 

na e Y^ g qps t j 


Jlps 


(1-e 2 ) li (Jl-2p+s)-(Jl-2p) Q 


jiqps ,f 


liqpsj 


4TTGp 0 


• ± k a * ♦! ' 


na ( 1 -e 2 )' 2 sin i Y pqps j 


2 Jl ♦ | / $P S 


[(A- 2p) cos i-q] Q JqP sJ 


...( 156 ) 


(1.57) 


(1 58) 


(1 59) 
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. (1.60) 


AQ 


sqpsj 


= 4TT G P° G ^ps QjlqpsJ / R x*+2/ 1*kj x 

1 a } V oo ♦ , / 


na ( 1 -e 2 )' 2 sin i 


Y pqps J 


(^)-f > (5-) 2 (s cos , a-2p)-4 cos 

j i ( 1 -e z ) y H 


oqps ,f 


• (Jl-2p +s) - q ) 


Aco„ , - 4tT G P» Q »qps J fR\ ,,z / l * t i \ 

„ 2 ATf * , ^ 
na( I -e 2 ) 2 0 jqpsJ 


— (^iEL.) F„p - cot , (iftt) Gjj ps » ^ - C2 ° nFw y»" (*) S 

e V } K d\ } 9?s 2 Q 2 0 - e 2) 2 . ± ^ q ^ 

0 pqpsj 


•[3 cos 2 i(Jl-2p)-q cos i + -f (1-5 cos 2 i)( 1 -e 2 P (Jl-2p +s)- -| (Jl-2p)j 


( 1.6 1 ) 


AM . _ 4TJ ^Po ^iqp ^pqps ,j ( R ( 1 \ 

' qps ~ ~p : ( ^ } ( w 

na 0 jiqpsj 


1 ^ ^ 

e 3e 


♦ 2 «*.,G ^ 9 CgpnGa , ( R .) 2 fs C0S 2 , w . 2p) 


^ P 5 '* (1_e2> Y^qpsJ 


- 2q cos i -U-2PM1-3 cos 2 i)( 1 -e 2 P (Jl-2p +s)] 


(I 62) 
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Perturbation amplitudes given by (1.57) (1 62) can also be 
expressed in the form of position perturbations in the radial 
(Ar), transverse (At), and normal (An) directions (Williamson 
and Douglas, 1972) 


i 2 

At = 9(1 + 5 e 2 ) V 2 (am + Alo ♦ AO cos 1 + (Ae 2 + (eAM) ) 2 ) 


Ar = \ Aa 2 + k 


(eAa + aAe) 2 + (aeAM) 


■ir 




An = a( 1 + ^ e 2 ) V 2 ( \ ai 2 + AO 2 sin 2 i)' 2 


(1.63) 
(1 64) 
(1 65) 


and for the velocity perturbations 


AV = 


/ Aa_ 's 2 


Ae _ Aae n 2 + ( ^ aM ) 
^ /a 2 a 3/ 2 ; 


/a 


, -> -v, dO 

a( 1 ^ e 2 ) 2 cos i 


2 ] / 2 


AY = 4 


Ae eAa f e ' 2 


/ ac 

^ /a 


2&> ' v /a 


am)' 


1 , 

"2 


(1 66 ) 


.( 1 67) 


and. 


AV = a( 1 + i e 2 ) 2 ' 

n 

/ 

\ 

2 

((-f 

' dt } 

2 A ,2 

+ n Ai 

+ (™) 2 sin 2 l • 


\ 

- 


/ 


(1 68 ) 


For formulas ( 1 63) through ( 1 68), we note the perturbation amplitudes 
are the RSS values of the perturbations and that small eccentricity 
approximations have been extensively applied 
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Formulas for the perturbations due to the solid earth tides can be 
obtained from ( 1 52) ( 1 57) by making the following substitutions 


OCEAN TIDES 


Replace 

Gp. R ( 
r -' ' 29 + 1 



SOLID EARTH 


With 

^ 2 J ( 


3- m \ 
3 > 


a 

q 




(1 69) 


From formulas ( 1 57) through ( 1 62) one can derive principal 
long-period perturbations for ocean or solid earth tides when 

P = ^2 - s = 0, + sign and q = m 


Then the angular argument becomes 


TlmVf = (i) 


(2-2h) co* ♦ (2-2h+_j)M* + kO* + mO + 71 


Tf p + 

^ O - + £ n L 

2 HmJ 


(1 70) 


with frequencies given by 


m%o,f 


(?) 


(2-2h)co* ♦ (2-2h+j)M* + kO* ♦ m6 


(171) 


It is apparent from (171) that the principal perturbations due to the 
solid earth and ocean tides have the same periods and therefore 
cannot be distinguished The satellite senses the total tidal effect. 
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APPENDIX B 


MATHEMATICAL FORMULATION 
OF THE SECULAR CHANGES 
IN THE LUNAR AND SOLAR ORBITS 


Due to t Me anelastic behavior of the Earth's response and the 
dissipation of energy in the oceans, the earth's tidal bulge 
causes secular changes in the motion of the moon and the sun 
The moon's orbit and the sun's apparent orbit can be treated in a 
way similar to that described in Appendix 1 and analytical per- 
turbations in their orbital elements can be obtained from equa- 
tions (1 57) through (1 62) it is evident, however, that secular 
perturbations will be generated only when the angular arguments 
given by equations (1 4) and (15) are independent of time 
Necessary conditions for the cancellation of the angular elements 
of the sun and the moon are that the motions of both the disturb- 
ing and the perturbed bodies, which are the same in this case, be 
referred to the same plane The mathematical formulation of the 
tidal potential given in Appendix 1, is in a system where the dis- 
turbing bodies are referred to the ecliptic whereas the perturbed 
body is referred to the equator The formulation of the potential, 
where the ecliptic is the reference frame common to both bodies 
but the tide amplitude and phase are in the usual equatorial 
reference frame, is given herein The formulation of the potential 
is developed for the case of the oceans only since the solid earth 
is a special case of this development. 


From equation (1 5), we can write: 

°SqJ r ^ ± qX 


where 


^ = (+)^(2-2h) oo*+ (2-2h+j) M* + kO*J 


+ m0g _ -I 


(2 1 ) 


(2 2 ) 
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Letting 


/ i + \ + 

A «,f = 4 nGR P. (^ 77 ) c J».f 


and noting that 


cos (<$ (±) qX ) = cos ( qX(±)^) 


the potential given by (13) can be rewritten as: 

U W,f = A Jq,f (rf" P !1 (S ' n 4) C0S 

which is the real part of the complex potential 

u *.f = A «.f (f) p *< s, "‘ )e 


where 


i = v-T 


Substituting the right ascension (cx) for the longitude ( X ) in 
equation ( 2 . 6 ), gives 


u pq,f = 6 


((£)*♦> Pw(sin#)e ",o<) 


Using the approach of Goad (1977), 


( R v Q + 1 _ . . * jqo< / R 7 P+1 — y y w<±> 

(-J P, q (sin ♦) e ' = (tJ (*-q>! v 4„ ¥ 


pq v 


• P^ v (sin u*) e 


v = o {±} 

[{±>vo<* - (v-q) %] 
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(2.3) 


(2 4) 


(2.5) 


(2 6 ) 


(27) 


(2 8 ) 



,<±> 


where the starred elements refer to the ecliptic and the function 
is given in equation (19) of Appendix 1 . Substituting (2 8) to (2 7) gives 

± t l | ^ e i(-«V±^ +(q ' v) T> 

*q,f - "j>q A (i-q)l v % g } 


YjqV ((-r/ +1 p pv(sin •*)? 


l({±}vo<*) 


(2 9) 


From Kaula ( 1 966): 


(frp )y (sin««)e' y “*=(f.ri I F )vp (i*) G !ps (e*) e 1 c J“ps 

p=0 S=-00 


..(2 10 ) 


where 


l* = ( j) - 2p)co* ♦ ( j)-2p + s) M* + vQ* - j mod ( U - v,2) (2 11) 


Ovps 


and therefore (2 9) becomes 


u pqJ = A i»q^ (i-q)i 


^ e i(-qOg(±) ^ + (q-v) t ^ ^ R. 


v>P,s 


F,v P 0*> G Jps (e*){Y * v e lC ?vps + Y 


jqv 


-i t* 

e pvps 


) 


( 2 . 12 ) 


or 


I j + x + 

U pq,f = A pq-f (P-q)l v,pX{ ± } K 8 


^ (4,/*' |: 1 .vpO*)G, ps (e«) 


{±} „l(< ± > C jvps" <,0 g (±) ^ +(q " v) ^' ) 


• ¥ e 

T jjqv c 


(2 13) 


This is the desired form of the potential, with both the disturbing 
body and perturbations referred to the ecliptic Note that the 
Apq f still refers to the usual ocean tide coefficients 
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The angular argument in equation (2 13) is given by 


r =(tU* vps - qG s (!) -3 * (q-v) y 


(2 14) 


or, expanding. 


r = {± } [( Q - 2p)oo* + ( D-2p ♦ s) M* + vO* - y mod ( j) -v,2)J -q0g ♦ (q-v) j 


+ [( + ) ((2-2h)co* + (2-2h+j) M* + kQ* ) + m0 g 


mTT 


+ 71 + 8 


j 


(2 15) 


The conditions tor secular effects are attained when the 
time-dependent angular variables in (2 15) drop out of the 
equation These conditions are 

1) -q0 a + rn0£ = 0, m,q > 0 ==> q = m + is + 

2) { + } vQ* + (?) kQ* - 0 ; v,k > 0 =» v = k, {+} = (+) 


3) { + } (jl - 2p) + (?) (2 - 2h) = 0 — * # = 2, p = h 

4) { + } ( j) - 2p + s) ± (?) (2 - 2h + j) = 0 ==» S-J (2 16) 


Note that when m = 0 or k = o, both sign alternatives exist, and 
we have applied the selection rule appropriate to the general case 
These degeneracies or symetries are the same as those discussed 
in Appendix 1 for (1.24) through ( 1 27). 

The only free variables for the secular effects are 
(±), h, j, k, m 
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wmch are also trie variables which define the Doodson argument 
in ( l 23) This implies that each Doodson frequency has a secular 
effect on the orbit of the sun or the moon due to the ♦ tides Any 
secular effect from the - tides are strictly through the m = 0 or 
k _ 0 degeneracies noted above, which, in the Doodson development, 
are incorporated into the Doodson coefficient amplitudes 


The constant phase angle r«*c associated with these secular 
effects is 


!>C = 1 ♦ 


~ mod 


k, 2) - 


krf 

o 


+ n + 


2m , f 


(2 17) 


which for 


k ^ 0 


for (+) 

V - \ | 

for (-) 

k r 2 


f%ec - n + E 2m j 


fsec - £ 


2m, f 


r sec = T1 "*■ £ 


2m, } 


Tsec - £ 


+ 

2m J 


(2 18) 


The n merely indicates a sign change and therefore the r S ec is 
nothing other than measure of the phase lag £ 2 + m j of the ocean 
tide potential 

We can also rewrite the amplitude of the potential in equation 
(2.13) as a function of the Doodson coefficient, A using equation 
( 1 7) of Appendix 1 


U 


♦ 


2m J 



moon 


sun 



(-o' 1 / 3-m (2*ro)i 

( 2- 6 0 m; 4 


e 


J Tsec 


(2 19) 
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The appropriate real form of the disturbing potential for the 
secular effects is the real part of the complex potential in (2.19) 
subject to the conditions of (2 15) Also, as equation (18) gives 


d = 6 0 k 



-k {±} k 


it can easily be shown that for any k, 

(- 1 ) d COS Tsec = cos £ ~ < 

2m, f 

Therefore the secular potential is expressed as 


(2 20 ) 


U 


+ 

2m J 



COS 



(2 21 ) 


where 


D 


2m, f 



moon 


sun 



( 2 + m) ! 
( 2 - m) 1 


( 2 . 22 ) 


Of particular interest are the secular changes in the semi-major 
axis, the eccentricity and the inclination of the orbits of the 
moon and sun These secular changes describe the evolution of 
the orbits of the sun and the moon and furthermore define the 
secular deceleration in the rotation of the earth, through the 
balance of the angular momentum equation 
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Given a potential U, the secular rates in the Kepler elements 
of the sun and the moon can be computed using the Lagrange 
Planetary Equations from Kaula ( 1966, equation 3 38): 


da* 

dt 

de* 

dt 


= S'nC.f 


2 ht 

<?) n." a *2e* l( ' ' f* 2 )’' 2 (2-2h*j) - (2-2h)] f sir, E 2 * m( 


di* 

dt 


= ( i) 0 ;„ , s,n E * 


n*a* 2 ( 1 -e* 2 )' 2 sin r 


f ~ t oil C ^ 

2m, f 2m, | 


(2 23) 


where n* is derived from Kepler s law: 


n* 2 a *3 _ u + 


When it is desirable to use the rate of change of the mean motion, 
dn* , m place of da* the appropriate relationship is 
dt dt 


dn* _ - 3 _ _n* der* 
dt ’ 2 a* dt 


The secular perturbations due to the solid tide are easily computed 
from equations (2 2 1 ) through (2 23) if the correspondences of ( 1 69) 
are applied 
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